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This article describes the frequency- acquisition performance of the Costas 
cross-over loop which is used in the Advanced Receiver II (ARX II) to perform 
Quadrature Phase Shift Keying (QPSK) carrier tracking/ The performance of the 
Costas cross-over loop is compared to two other QPSK carrier tracking loops . the 
MAP estimation loop and the generalized Costas loop. Acquisition times and prob- 
abilities of acquisition as functions of both loop signal-to-noise ratio and frequency- 
offset to loop-bandwidth ratio are obtained using computer simulations for both 
type-II and type-III loops. It is shown that even though the MAP loop results in 
the smallest squaring loss for all signal-to-noise ratios, the MAP loop is sometimes 
outperformed by the other two loops in terms of acquisition time and probability. 


I. Introduction 

Quadrature Phase Shift Keying (QPSK) is an efficient 
modulation scheme which has been used extensively in 
communications systems. It offers a significant advantage 
over Binary Phase Shift Keying (BPSK) since it permits 
the transmission of twice the data rate in the same chan- 
nel bandwidth [1], as both the in-phase and the quadrature 
components of the carrier are modulated. An even better 
use of bandwidth can be obtained by going to higher-order 
phase shift keying, but that would require a higher signal- 
to-noise ratio to achieve a required bit-error probability. 


Depending on the application, derivatives of QPSK can 
also be used; these include Staggered QPSK (SQPSK), 
which is better suited for satellite channels with Travelling 
Wave Tube (TWT) amplifiers and Minimum Shift Keying 
(MSK), which can be thought of as SQPSK with sinusoidal 
pulse shapes. 

One of the functional requirements of the Advanced 
Receiver II (ARX II) is that it acquire and track QPSK 
signals. The most likely candidate algorithm for the ARX 
II is the Costas cross-over loop because of its ease of imple- 
mentation and relative performance. Carrier synchroniza- 
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tion for QPSK has been addressed by many authors in the 
literature and several loops have been introduced and an- 
alyzed in terms of their respective phase-error variance [2- 
4]. However, very little has been reported on the acquisi- 
tion performance of these loops, mainly because the prob- 
lem becomes highly nonlinear and precludes any analysis 
that provides insight into loop behavior. The suppressed 
carrier tracking loops are typically derived from Maximum 
A Posteriori (MAP) estimation theory; the derivative of 
the likelihood function is used as an error signal to pro- 
vide a closed-loop implementation commonly referred to 
as MAP estimation loops [2]. 

Three QPSK loops are considered in this article in 
terms of their acquisition performance: the MAP esti- 
mation loop, the Costas cross-over loop (also known as 
the polarity-type Costas loop), and the generalized Costas 
loop (sometimes referred to as the conventional Costas 
loop). Mean acquisition times and probabilities of acquisi- 
tion as functions of both loop signal-to-noise ratio (SNR) 
and frequency-offset-to-loop-bandwidth ratio are investi- 
gated using computer simulations for type-II and type-TII 
loops. 

The loops are all shown in Fig. 1 where different pro- 
cessing of the in-phase and quadrature channels identifies 
which loop is being implemented. The MAP estimation 
loop produces an error signal given by 

Z 0 = Z c Unh(RZ $ ) - ZMnh(RZ c ) (1) 

where R is the symbol SNR and Z C) Z s are the in-phase 
and quadrature signals respectively. This loop is “diffi- 
cult” to build because it requires a priori knowledge of the 
symbol SNR and implementation of the hyperbolic tan- 
gent operations. The latter is hard to achieve in an analog 
fashion, but is easily done digitally with current technol- 
ogy using read-only memory (ROM) chips. To reduce the 
complexity of the MAP estimation loop, approximations 
to the hyperbolic tangent function can be used at low and 
high symbol SNRs. For instance, at high symbol SNRs, 
the hyperbolic tangent can be approximated by the signum 
function (tanh(x) ~ sgn (x)) to yield the Costas cross-over 
loop with error signal 

Z 0 = Z c sgn(Z,) - Z,sgn(Z c ) (2) 

At low symbol SNRs, the nonlinearity hyperbolic tangent 
can be expanded using a Taylor series to give tanh(x) ~ 
x — x 3 /3, which results in the generalized Costas loop with 
Z 0 given by 

Zo = Z C Z S (Z 2 - Zl) (3) 


The two approximate loops are straightforward to im- 
plement and do not require knowledge of the symbol SNR. 
However, all three loops require an estimate of the incom- 
ing signal amplitude in order to operate the loop at the 
desired bandwidth. 

Section II of this article develops a mathematical 
model for the various signals in Fig. L Simulation results 
are presented and discussed in Section III followed by the 
conclusion in Section IV. 

II. Mathematical Model 

In a typical communication system, the received sig- 
nal is first downconverted to an appropriate intermediate 
frequency (IF) for further processing. At that point, the 
QPSK signal can be modeled as 

r(i) = y/S[d](t)s\n(uJit 4- 9) + dq(t)cos(u>it + 0)] + n (0 

(4) 

where S is the received power (watts), u the IF radian 
frequency (rad/sec), 9 the signal phase (rad), and dj(t) ) 
dq(t) the in-phase and quadrature data streams given by 

+ oo 

di(t)= E a Ki p{t-kT) i = I,Q (5) 

fc= — oo 

where p(t) is the baseband Non-Return-To-Zero (NRZ) 
pulse limited to T seconds and a kt i the equally likely ±1 
binary symbols. For SQPSK, the in-phase and quadrature 
baseband pulses will be offset by half a symbol period and 
for MSK, the pulses will be sinusoidal. The narrow-band 
noise process n(f) can be expressed as 

n(t) = \/2n c (t)cos(uJit + 0) — y/2n a (t)s\n(uit + 9) (6) 

where n c (t) and n,(£) are statistically independent sta- 
tionary band-limited white Gaussian noise processes with 
one-sided spectral density No (watts/Hz) and one-sided 
bandwidth W (Hz). Assuming a digital implementation, 
the signal is first digitized at a high rate of R s samples per 
sec ( T s = 1 /R s is the sampling time) and subsequently 
demodulated using the in-phase and quadrature references 
(see Fig. 1). The resulting samples are then accumulated 
over a symbol duration and normalized to produce 

Z 9 (iT) = aijco$<f> — ai t qsm<j> — Njcos<f> — Nqs'iiuP (7) 

Z c (iT ) = aijsin<f> + a^gcos^ — Nj s\n<f> + Nqc os<p (8) 


151 



where <p - 6 - 0 is the phase error and Nj, Nq are zero- 
mean Gaussian random variables with variances Nq/TS. 
It is assumed that perfect symbol synchronization has been 
achieved and no attempt is made to quantify the results 
of timing jitter on the carrier-acquisition process in this 
preliminary study. However, the acquisition performance 
with unknown symbol epoch is of prime interest as it re- 
flects the carrier-acquisition process in a practical system. 
Moreover, the various loops might react differently in the 
absence of symbol synchronization and a complete study 
incorporating those effects is necessary. Depending on the 
loop, the error signal Zq is computed from Z c and Z 3 using 
the appropriate equations as discussed earlier. Further av- 
eraging over N symbols can be performed before updating 
the loop filter, which in turn controls the frequency of the 
numerically controlled oscillator (NCO). 

The tracking performances of these loops have been 
derived elsewhere and are repeated below for convenience. 
In all cases, the variance of the phase-error process is given 
by 


where p is the loop SNR of the QPSK loop, p c = S/NqBi 
the loop SNR of a “classical” phase-locked loop [1], and 
Si is the “squaring” loss which for the MAP estimation 
loop is given by [2] as 


loop outperforms the other two loops for all values of R. 
For R > 8 dB, the performance of the Costas cross-over 
loop approaches that of the MAP loop. Similar arguments 
apply for the generalized loop at low SNRs (R < —4 dB). 
The analytical results were reproduced through computer 
simulations to obtain greater confidence in the programs 
before generating new results on acquisition. 

III. Simulation Model, Results, and 
Discussion 

The loops in Fig. 1 are simulated by generating the in- 
phase and quadrature signals Z s and Z c according to Eqs. 
(7) and (8). Depending on the loop under consideration, 
the instantaneous error signal Z 0 is computed using Eqs. 
(1), (2), or (3). The simulations were carried out assuming 
N = I, i.e., the loop update rate is identical to the symbol 
rate. In this case, the NCO transfer function is given by 

N(z) = ( 13 ) 

(where T u denotes the loop update time) and it relates 
the filter output to the incoming phase estimate in the Z- 
domain. In some digital implementations, the NCO might 
have an additional delay, in which case the loop perfor- 
mance will be slightly worse than the one considered here. 
Unlike tracking, which depends solely on the loop band- 
width, acquisition behavior is a function of the loop filter 
itself and not just the loop bandwidth. The simulations 
were carried out with the following filter 


[tanh(fl - y/Rx) - Rsed^jR - v/fix)] 2 
(1 + #)tanh 2 (ft - V Rx ) - [(x - \fR)tanh(R - VRx)] 2 

( 10 ) 


where x is a zero-mean, unit-variance Gaussian ran- 
dom variable and the overbar denotes statistical ex- 
pectation. For the Costas cross-over loop [3] 
(erf(*)=(2/v / i r ) Jq e'^dt denotes the error function), Sl 
becomes 


q _ [ e rf(y^72)-y^RA ^ B/ T nn 

^ L l + fl-^2/* e-*/ 3 + '/7? erf(,/;?/2)P ^ ' 


and 


5 ^ 


1 

1 + 2k + 7 ?* + 2^ 


( 12 ) 


for the conventional Costas loop [4]. In all cases, Si is 
a function of only one parameter, the symbol SNR R de- 
fined as ST /N 0 . Figure 2 depicts the “squaring” losses 
for the various loops; it is clear that the MAP estimation 


F{z) = G 1 + t%=t + (14) 

where 

Gi = rd/T u (15) 

G 2 = rd 2 /T u (16) 

G 3 = krd 3 /T u (17) 

and 

d = 4B L T u (r - Jb)/r(r - k + 1) (18) 

Bi denotes the design-loop bandwidth in Hz, r is typically 
2 or 4 and is equal to 4£ where £ is the damping ratio, and 
it is a type-III loop gain parameter (k = 0 for a type-II 
loop) with typical values ranging from 1/4 to 1/2. The 
filter of Eq. (14) was derived from an equivalent analog 
filter using the Impulse Invariant Transformation (IIT) [5]. 
The actual loop bandwidth might be larger than the de- 
signed Bi depending on the product BiT u . Generally, for 
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BlT u < 0.05, the actual loop bandwidth is very close to 
the design loop bandwidth B L . 

Each simulation point was obtained using 5000 runs 
for each combination of loop SNR and frequency offset. 
Each run is for a maximum of 100 [Bl seconds and the 
loop is declared to have acquired the incoming signal if the 
magnitude of the residual phase error becomes less than 
tt/ 8 radians and remains there for at least 10 [Bl seconds. 
By residual phase error is meant the phase error modulo 
7 t /2 since the stable lock points are integer multiples of 
7r/2. The “acquisition” definition used in this article might 
be too strict for some applications, especially at low loop 
SNRs. This will be discussed later in more detail. 

Since the acquisition results are based solely on com- 
puter simulations, it is worthwhile at this point to discuss 
their confidence level. Note that the outcome of each run 
is either “loop acquired” or “loop did not acquire,” i.e., 
a binary decision. Let X be a random variable denoting 
the number of “loop acquired” decisions in n runs. Then, 
X can be modeled by a binomial distribution with mean 
np and variance np(l - p) where p denotes the expected 
probability of acquisition in a single run. Hence, the error 
level (ratio of the standard deviation to the mean) becomes 
[(1 — p)/np] 0 5 (accuracy level is 1 — [(1 — p)/ n p] 0 5 )* As 
an example, when simulating the loops with an expected 
probability of acquisition of 10 percent using 5000 runs, 
the results will be correct with a 95. 4-percent accuracy 
level. On the other hand, for a 90-percent probability of 
acquisition, the accuracy increases to 99.53 percent. 

Probability of acquisition as a function of time (nor- 
malized by loop bandwidth) is depicted in Fig. 3 for all 
type-III loops and in Fig. 4 for all type II loops. It is 
clear that for a given acquisition time, the acquisition 
probability increases with decreasing frequency-offset-to- 
loop-bandwidth ratio A //Bl, as expected. Furthermore, 
type-II loops seem to outperform their respective type-III 
loops in terms of acquisition time for AJ/Bl < 1, but 
type-III loops are comparable for A //Bl > 1 for “high” 
acquisition probabilities. As an example, a type-II gen- 
eralized Costas loop requires about 18/ Bl seconds to ac- 
quire a 0.5 B L -llz frequency offset with a 60-percent prob- 
ability, whereas the corresponding type-III loop requires 
about 30 /Bl seconds for the same probability. On the 
other hand, when the frequency offset is about 1.25i?£,, 
then the type-II and type-III loops require about 90/ B L 
and 40 /Bl seconds respectively for a 20-percent acquisi- 
tion probability. This is better shown in Fig. 5 for the 
type-II and type-III Costas cross-over loops where either 
type can dominate depending on the ratio A f /B l* Note 
that at high probabilities of acquisition (> 90 percent), the 


performance of type-II and type-III loops will become com- 
parable as shown by the trends of Fig. 5 for A// B L — 0.75 
and 1.25, with type-II loops seemingly in the lead. 

When comparing the different loop structures, the 
MAP estimation loop is sometimes outperformed by one 
or both of the other two loops; for example, at 50 percent 
acquisition probability with a frequency offset of 1.0 Bl 
Hz, the MAP estimation loop and the Costas cross-over 
loop require about 73/ Bl and 68/ Bl seconds, respec- 
tively. This is an unexpected result but not very sur- 
prising as the MAP loop is the optimum closed-loop es- 
timator for an unknown phase but not necessarily for an 
unknown frequency. Figures 6 and 7 compare the acqui- 
sition performance of the three different loops structures 
at IS- and 25-dB loop SNRs for both type-II and type- 
III loops. It seems that for type-II loops, the MAP loop 
is outperformed by both the Costas cross-over loop and 
the generalized Costas loop, depending on the loop SNR 
(Fig. 6); however, for type-III loops, the performance of 
the MAP loop is almost identical to the best of the other 
two loops (Fig. 7) for all practical purposes. 

Another parameter of interest is the probability of 
acquisition assuming an unlimited acquisition time (sim- 
ulations were actually carried out with 100 /Bl seconds). 
This differs from the cumulative probability of acquisition 
plotted earlier which employs a fixed acquisition time. The 
probability of acquisition as a function of frequency-ofiset- 
to-loop-bandwidth ratio for a given loop SNR is plotted 
in Fig. 8 for type-II loops and in Fig. 9 for type-III loops. 
As expected, the probability increases with increasing loop 
SNR but decreases with increasing A f /B l- It is interest- 
ing to note that for both type-II and type-III loops, the 
MAP estimation loop for all practical purposes achieves 
the highest probability of acquisition for a given loop SNR 
and a fixed A//J5 l- However, from the previous discus- 
sion, the MAP loop does not necessarily have the highest 
probability of acquisition for a fixed acquisition time. 

Finally, Fig. 10 depicts the effect of B l T u on the cu- 
mulative probability for the Costas cross-over loop. As 
expected, better performance is obtained by decreasing 
BlT U) which results in a faster loop update rate for a 
given loop bandwidth. 

As mentioned earlier, the acquisition definition in- 
voked so far might be too strict at low loop SNRs. This 
is because even though the loop can acquire the incoming 
phase (i.e., the residual phase error becomes less than tt/8 
as required), it cannot remain there for 10 /Bl seconds be- 
cause the mean time to lose lock at low loop SNR becomes 
“small” and affects the probability of acquisition. This is 
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illustrated in Fig. 11, where the cumulative probability of 
acquisition is plotted against Bit for different lock criteria 
for the Costas cross-over loop at both 15- and 20-dB loop 
SNRs. For a 20-dB loop SNR, the cumulative probability 
of acquisition is for all practical purposes independent of 
the lock criterion used because the mean time to lose lock 
is larger than the 10 jB L seconds required. However, for a 
15-dB loop SNR, the mean time to lose lock is small com- 
pared to 10 /Bi seconds and the performance varies with 
each lock criterion as expected. In that sense, the lock 


criterion of 10 /Bl seconds might not be fair for low loop 
SNRs. 

IV. Conclusion 

Three QPSK carrier tracking loops were compared in 
this article in terms of their acquisition performance. It 
was determined that the Costas cross-over loop is the most 
likely candidate for implementation in the ARX II due to 
its simplicity and overall relative performance. 
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CUMULATIVE PROBABILITY OF ACQUISITION 



Fig. 3. Cumulative probability of acquisition versus B L t for type-ll 
loops: (a) Costas cross-over loop, (b) generalized Costas loop, 
and (c) MAP estimation loop. 


CUMULATIVE PROBABILITY OF ACQUISITION 



Fig. 4. Cumulative probability of acquisition versus B L t for type- 
ill loops: (a) Costas cross-over loop, (b) generalized Costas loop, 
and (c) MAP estimation loop. 








Fig. 5. Cumulative probability of acquisition versus B L t for the 
type-11 and -III Costas cross-over loops. 




B L t 


Fig. 7. Cumulative probability of acquisition versus B L t for all 
three type-ill loops: (a) p = 15 dB and (b) p = 25 dB, 


Fig. 6. Cumulative probability of acquisition versus B L t for all 
three type-ll loops: (a) p = 15 dB and (b) p = 25 dB. 
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PROBABILITY OF ACQUISITION 



Af/B^ 


Af/B^ 


Fig. 8. Probability of acquisition versus A f/ B L for type-!! loops: Fig. 9. Probability of acquisition versus A f/B L for type-ill loops: 

(a) Costas cross-over loop, (b) generalized Costas loop, and (a) Costas cross-over loop, (b) generalized Costas loop, and 

(c) MAP estimation loop. (c) MAP estimation loop. 
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